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Abstract 


(DRIRU) gyroscopes on the Upper AtLsphere Research Sat nt niAR^" 6 ^ r °*° r ' nertial reference unit 
(EUVE). The accuracy of „ t'T 1 a " d th ® Extreme U,traviolet Explorer 

spacecraft attitude is sensitive to gyro noise EUVE nvm Ltl ^®° nthms that , use W° data for Propagating the 
model, which is used onboard vahous spacecraft VaHdate 3 sin 9 le - axis ^ ro noise 

effects on the gyro noise model and the overall impact on^ttitnH^nT 568 ^ potential Impact of temperature 
density (PSD) of the gyro noise is m flioM data ISTTF , The P ° Wer spectral 

of actuator dynamics on the PSD function is also discussed. 9 * by F 1 Four,er Transf °rm (FFT) The role 

Introduction 

c^he^ro ^o^ 0 This over a period of time are affected by the 

Atmosphere ReseSch Satellite OJAR^an^pi Chara * enze that noise using in-flight data for the Upper 
have Ledyn^otf ^ ^ 

specifications of this hardware are presented below ^ ^ gyrOSCOpes on 5031(1 The description and 

^£?js modei that is — *■ ^ 

EUVE and I JAR ^ fnr T ^ board vanous spacecraft (the Gamma Ray Observatory (GRO) 

m™re d „ts U o^i ™ ^ Sl "“ EWE “ 

“ rs " o ;r a ,™; t T 0 ,s ^ y c : c *r ra “ 

- 

instrument operation. The power spectral densitv rpcnwth easurement noise may be due to science 
Transform (FTT) to evaluate UOs suasion Z' m ® ” 0ISe •* obtained b > Fas ' four.cr 

•ha, iso„, y me ,L noise «*" — » * *« 
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IRU Description and Specifications 


The DRIRU II consists of three gyroscopes, each ^ a total of six channels of 

degrees of freedom and rate information along two body _ ( temperature of each is reported in the 

information. Each gyroscope has internal temperature each body axis. When the 

downlinked telemetry. The gyroscopes ^e confirm f eorie nted to maintain a null deflection from an inertial 

spacecraft attitude is noninertial, the IRU gun t0 accoin piish this reorientation is proportional to 

orientation. The current required to produce current is converted to a series of pulses, which are 

the angular rate about the corresponding axis. q provides rate information as incremental 

countefand reported as accumulated rotation So accumulated angles reported every 

angles every 64 milliseconds and a separate electro every 128 milliseconds and 

- — - - — 

together as the IRU. 

operate two rate rap f a The 

low-rate mode allows rates of up to 400 arcsec/see <T1 ^ „ fo , the 

the ordro of a, leas, 6 hours,. Only dte short-tenn noi* perfonnauce ,s 
pertinent for the onboard attitude determination processing. 

The mil accumulated rotation* angle meatmen, safe reported SS S 

as. rr », ^ - 

factors, and biases. 


Noise Model and Rate Noise 


A common model for IRU rate noise (documented in Reference 2 and used in the 
attitude determination) uses the sum of two noise processes: 


onboard Kalman filter for 


r(r) = v(r) + MO 


, f ,• „ ranidlv varying random process that is modeled by ® 

where r(r) is the ra ^°^ “ * slowtyvar^ng random process that is modeled by a random walk process (see 
white-noise source, and b(r) is a s y ary g fornue noise and the bias white noise, «(r), 

“dm g t^^ 

iidjua^for'hi^^ p^« ^(Tjcan be taketTas zero a, dm «. Ume. % . of the timerpan under 

consideration The autocorrelation of the white-noise source has the following form. 

(v(r)v(r')) = <^r -r') 


where a] has the dimensions of angle squared per time 
random walk process is the integral of another independent 
that drives the random walk of the rate such that 


and <...> is the statistical expectation operator. The 
white-noise source. Let u(t) be the white noise source 


(u(t)u(t')) = crls(t-t') 
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per time cubed. Then, the random-walk process is then given by 

t' = t 


where cr has the dimensions of angle squared 

Ht) = r 1 u{t')dt' 

Jf = t 0 


(b(t)b(t')) = oi min((r-/ 0 ),(r'-r 0 )) 

where min((r-r 0 ),(r'- r 0 )) is the smaller of t-t, 
autocorrelation of the rate noise is given by 


o and t'-t 0 . For independent random processes, the 


(r(r)r(r')) - <f6(t-t') + </min((r-r 0 ),(r' -r o )) 

Observable Consequences of Rate Noise 

The rate output of the IRU for each channel has the following form 


n(r) = <y(r) + r {t ) 


where fi(t ) and o»(r)are the observed and 
angle 


true rates, respectively. The IRU converts this 


rate to an accumulated 


©(/) - J L2(r) dt = J <y(r) di + j r(t ) dt 

Let the rate be integrated from time r 0 . Let q(t) be the noise of the integrated rate as a function of time Then, 

q(t) = si = i 0 r{t ' )dt ' = si ~i v(t ' )dt> +si *(»')*' 

The variance of the noise of the angle accumulated from time t 0 is 

ft' = t 


( q{t)2 ) = Sr = l dt Sr = l dt ‘i( v (t') v ('")) + (b(t')b(t"))) 

7 / x rt f = t c t u = t 

^ ' to ' + Jr = t 0 df 'J r = t / r ^ min(t '- t O’ r - t o) 

= ° 2 J'-t 0 ) + j° 2 Jr-t o y 


(1) 


^^ f,Ca .K°7 fth , C 1 b °i Ve equatl0n is the 8° al of our gyo noise processing. For EUVE all six gyro channels are 
reponed ,hc downlinked «le m e,ry Th, s luxury allows us observe rhe £ro no.se direct If 
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me used where the spa«craft is undoing — 

alignments of the gyro channels are neglected, * h lcft ^ a quantity equal to the difference 

for a given axis will be to subtract the true rate ^Mon. On ms ^ then lett ^ a q ^ # ^ rf pure 

of the two channels' rate biases and rate noises. Ifthe ini (RMS) strength of 

noise is left that theoretically is described by the above equation, but with the root mean squa 

the random walk sources for the gyros given by 


( 2 ) 


and 


<t m /(V2) 


, , . . . c rSft-t ) 2 which is due to the 1-sigma uncertainty, 

An additional quadratic term can be included in Equation (1), cr b \t 0 ) > 

in the initial bias calculation. The bias uncertainty is discussed later in this analysis. 

Results from Applying the Noise Model to Observed Data 

Both low-rate mod. and high-mte mod. sym noise SESU « 

ISW «« ,L“t - 

difficult, eventually timespans when EUVE was m the survey phase were used exclusively. 

In Figure 1, the square of the propagated error due to the note 0 ^^ M^T^^Ws 

channels is shown as a function of of die variance plot over time (Figure 1) versus 

time, the gyros were in the Tow-rate mode Company ^ mght be by 

the temperature plot (for the primary Xchanne '^ ^ f the taken on that day, the apparent 

temperature changes. Further investigation into ^re^t to ffie Average (over the prime and backup 
bias (computed at 5-mmute interv Is), ^ r ° u ^ the averag e difference between two gyro signals for a 

channels) temperature as shown m Figure 3. As the temperature-variation induced effects From the 

given axis, both gyros must ^have constan P— rs, such as the gyro scale 
data processing performed for this analysis, it is unclear wheUier otn gy * m ^ strengths , are also 

factor for conversion of the telemetered counts to enginee 8 ExDlorer (COBE) Reference 3, has shown 
functions of temperature Previous analysis for the the DRJRU II. The 

gyro parameter dependence on gyro tempera res approximately 320 seconds after a change m gyro base 

COBE experience was that a variation in drift rate was PP s corre lated to the gyro temperature 

plate temperature Note that the bias .^^Jure^l is the integral with respect to time of the bias. Further 

fpta„.dTcC«. e ni , P h. relationship of >h. changes in ,he DRIRU 11 parameters wilh temperature 

variation. 

SSS £ SSTSw mlm ZSm over 

Sioni in time Fo, this reason, only EUVE some, mode dam were used. 

t , „ Arnn tc n? pfArpnee 4^ to fit the gyro noise model to on-oibit data, is 

rrriTf" 

the initial bias to the propagated error is analyzed in detail 
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Figure 1. X Propagation Error Squared for December 11, 1993 (Low Rate) 
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Figure 2. Prime and X and Y Channel Temperature for December 11, 1993 



Figure 3. X Channel Bias Versus Average (Prime and Backup) Temperatures 
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by the rate white noise, with strength denoted by <7 V . , 2 / 200 ("sec) This results in the 

variance is given by a* '300(sec), with the resultant bias vanance, a 2 , given by oj /300(sec). This 

final equation for the variance of the noise of the angle accumulated from time t 0 given by 


(q(t) 2 ) = o>[(t-t 0 ) + {t-t 0 ) 2 / 300(sec)] + 1 (r - r 0 )' 


(3) 


variance was assumed to be of the form 


Q(Ar) 2 = a, (At + Ar 2 / 300sec) + a 3 (Ar) 3 


(4) 


parameters over all the gyro channels. The final procedure is described below: 

, . Find limespans of data that are a. a near!, constant rate fo, Out entire Umespan; for EUVE, we rned the 
playback survey mode data from September 5 through September 11, 1992. 

? Fncure that there are no extreme temperature variations of the gyros. By visual inspection of the data, the 
iemnerature varietTno °niore ihan 1 Censor conn, (.14 deg Celsins) front its median value for the d« 
chosen (with the exception of infrequent, short-lived peaks of 2 counts, which are possibly noise 
temperature sensor output) 

3 Process ihe pnman. and backup channel information for each asds, differencing the backup channel from 
rp S r; Ze, ,0 remove the signal and leave o^y the value of the back^ channels no.se 
subtracted from the primary channel's noise and the difference ,n the two channels biases. 

4. Compute the initial bias using a 5-minute timespan and compensate for the initial bias throughoul the 
timespan of data used. 

The cubic equation describes the model for Ihe vanance of the 

the data taken can then be calculated and compared to expected values 

£22 SESSSSpSSSSS 

4. Data points are indicated by the squares; the fitting curve is the dark line. 

The 1-sigma white noise strengths for an individual channel that corresponds to the cubit fit shown above are 


<r„ = 0.12 arcsec/ sec 


1/2 


and 


(5) 


a u =5.21x10 5 arcsec/ sec 


1/2 
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Figure 4. Average Variance Over All Channels (High Rate) 

Typical values used for the DRIRU II in low-rate mode, (for GRO and EUVE, for example), 


are 


and 


< j v = 0. 042 arcsec/ sec 1/2 
< j u = 4.44 x 1(T 5 arcsec/ sec 372 


( 6 ) 


Our values are for the DRIRU II in high-rate mode, so no direct comparison can be made except that the computed 
values are higher and on the same order of magnitude This comparison to the onboard numbers which JKSd 

inlTiLTh for the '° w ; rate mode ’ IS made mainly to verify our results. Since our results are for the DRIRU 

, 8 6 ^ the EUVE and GR0 values are for the low-rate mode, attempting to draw any further 

rlT nS ab t t Ut the aCC “ ra 2 0f the onboard numbers would be specious because the parameters could differ 
KalmaiTfilter wh Basedon these data - the onboard white noise strengths used by the EUVE onboard 

Were m hl8h ' rate m0de Sh0U,d have ten mcreased from the ^ values^ 

Discussion of Results 

the e DRmu n’il°h Val h dat f ,h 17 Se r del USing on -° rt,it data ’ can be considered accomplished, although only for 
‘I™ hlg rate mode and for EUVE in surve y mode so that temperature effects are reduced The cubic 
onK^,ri aVerage C0V p ana , nCe 0f the propagat)on erTor due t0 noise is accurate enough from a practical sense for 
P^^ng Furthermore, early in the mission, when EUVE was using the gyros in the high-rate mode 

Lode the * urvey phase) > the gy f0 noise Parameters used by the onboard Kalman filter were for the low-rate 

mode. Based on the white no.se strengths computed for the high-rate mode, this could cause a decrease d The 
accuracy o t e onboard attitude determination algorithm, and was, in fact, seen (Reference 5). Based on 
recommendations from the FDF, the low-rate mode was used when possible for EUVE to improve the onboard 

the W ° U,d h3Ve ^ ,0 ^ ** n0iSC Paramders t0 the when 

•h^ Ver ’ ,he difficuIties due 10 chan ging gyro temperatures do cause concern. As discovered during our analysis 
3PPear t0 bC temperature dependent to a degree that might impact the usefulness of the noise 
model Further analysis is necessaiy, using low-rate data when EUVE is menially fixed, to quantify the impact of 
the temperature dependence on the noise model accuracy for the low-rate scenario. 
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Potential topics for further analyses are as follows. 

1 . Investigate the temperature effects on the gyro parameters 

2. Continue analysis of the low-rate mode to determine the impact of the temperature dependence on the 
applicability of the noise model 

3 . Process further data and compute the statistical confidences of the results 

Power Spectral Analysis 

7. 

For a continuous function, h(t), there exists a Fa f ¥o ^ f^ 0 n dSing 

SET ^..d limited. The Nyquis, frequency is defined ns (HR), where T is 

the sampling interval. 

If the function h(t) is continuous, then a sample of Jiff) at times separated by T can be represented as 


Y^h(t)S(t-nT) 


^ is o y f 

finite duration. The truncation (or window) function is defined below for a given ura ion 0 . 


*(r) = 1 -772 < t < T 0 -772 

= 0 otherwise 


( 8 ) 


Thus, the expression h s (t) 


h s (t) = ^ h(t)S(t - nT)x(t) 

H=-« ® 


( 9 ) 


represents N = TJT sample peints of the function W The tnmcalron of the time senes introduces rippling in the 
FFT. The discrete FFT of truncated sample wave form is given by 


H (n / NT) = J h{t)S(t - nT)x{t)e^"‘ 


( 10 ) 
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We use the FFT algorithm, popularized by Cooley and Tukey, to compute the FFT of the truncated sampled gyro 
signal. To obtain the power spectral density function of h s (t), suppose that the number of observations N = 2q + 1 
is odd. If we use a Fourier series to model the time series, such as. 


z(t)-a a + ^ {a, cos(2 >r/r) + # sin(2 nf t t)} + eft) (H) 

i=i 

where f i (= i/T) is the i lh harmonic of the fundamental frequency (1/T) and eft) is the random noise with zero 
mean and constant variance cr 2 The terms in the curly bracket represent the signals produced by the solar array 
movement and the operation of the science instruments. A least squares estimate of the coefficients a 0 , a J and /?, 
yields for i = 1 through q 


a 0 = gyro bias <z(T)> 

N 

a i = (2 / N) ^ z{nt)cosf2nf i t) (12) 

m - 1 

N 

b i = (2 / N) ^ z(nt)sin(2 ztfjt) 

m~ 1 


The spectral power or amplitude in the i 1 * 1 harmonic is 

P(fi) = (N/2)(a? +bf) t = 1, 2, 3,...,q (13) 

and result for even N is similar except that the b i are zero. The FFT of the gyro rates would yield the harmonic 
content of the gyro rates, i.e., the a v and, subsequently, the power spectral density. The expectation value of the 
variance of the time series is related to the PSD in a simple way. 

2 wo-{z(o)) 2 = X />(/;.) (i4) 

l i=l 

If the times series consisted only of white noise, then the amplitude of P{f) would be equal to 


PW-2a* (15) 

and its amplitude would be independent of frequency White noises ordinarily defined as noise that possess a flat 
power density spectrum for all frequencies. Evidently, if the PSD has a constant value for all frequencies, the total 
power represented by the noise would be infinite In practice, we defined white noise as a flat power spectral 
density over the frequency range of interest. However, any harmonic content in the times series would add to the 
value of variance at those harmonic frequencies. 

P(f,) = 2cf + N{a? +bf)/2 (16) 

In practice, it is not very' probable that the frequency' of the unknown sinusoidal component would match any of the 
calculated FFT frequencies. More likely, the variance would be spread among several frequencies and resemble a 
peak with finite width more than a delta function spike. 

The PSD is related to the autocorrelation R^t) by 


PSD= r R x (t)e 2 mfi dt 

J - oo 


(17) 
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We assume the autocorrelation to be described by 


(v(f)v(r’)) = cP~ v 


(18) 


and perform the integration to obtain 


PSD = a 2 


(19) 


where the units of c? v are angle 2 /time (arcsec 2 /sec). The PSD will have the units of the quantity being analyzed 
(rad/sec) squared per freq (sec’ 1 ). So, our units should be rad 2 /sec 2 . However, we are performing the discrete 
Fourier Transform. The relation between the discrete Fourier transform (H n ) of a set of numbers and their 
continuous Fourier Transform (H(f n )) when they are viewed as samples of a continuous function sampled at an 
interval dt can be rewritten as H(f n As we are computing the discrete Founer Transform of the rates 

and using the square of the magnitudes to represent the PSD, we need to multiply by dt 2 squared to get to an 
approximation to the continuous PSD, which is related to the continuous autocorrelation (with the delta function) 

as above. 


For low rates, the 0.042 arcsec/sec 1/2 value for the white noise and the equation 

P(f) = 2c 2 


gives a value for our discrete PSD of 


( 20 ) 


PSD = 2-crJ /(.128sec) 2 

= 2 (57.296) ' 2 -3600" 2 0.042 2 /(0.128) 2 (21) 

» 510~ 12 

Two criteria are used to identify significant peaks in the power spectral density. First, only peaks that have 
amplitude that are at least an order of magnitude larger than the background level are considered This criterion 
establishes the signal peak as statistically significant. Secondly, signal peaks that have finite bandwidths are 
indicative of complex physical processes 

The two main practical limitations for applying the power spectra method are as follows: (1) the length of the times 
series should not be a multiple of the frequency of interest and (2) the frequency of interest must be below half of 
the sampling frequency. However, in most cases, one does not know the frequency of interest. The most reasonable 
way to analyze the power spectra is to obtain spectra for different sample lengths. In this way, the occurrence of 
false periodicities are minimized. 

Power Spectra of UARS Gyroscopic Rates 

Data from UARS were used for two basic reasons: (1) the solar array motion was known to induce motion in the 
spacecraft and (2) some of the scientific instruments were suspected of causing motion in the spacecraft A search 
was made of the archival data for timespans in which investigators could isolate noise associated with the quiescent 
spacecraft from that of the operation of the solar arrays, the science instruments, and nominal operations For the 
quiescent period, June 4, 1992, was selected since the science instruments were turned off, and the solar array was 
parked On August 8, 1993, there was an extended period when only the solar array was in operation. Three days 
later, on August 11, at 17:50 Greenwich mean time (GMT), the UARS science instruments were turned on. 
February 6, 1994, was selected to represent nominal operations. 

On board the UARS spacecraft, the gyroscopes' digital rates are sampled every 0.128 sec, which is subsequently 
sent to the FDF for ground processing. On the ground, the digital rates are converted to gyro rates in units of 
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rad/sec. These data are ordinarily used to determine the spacecraft's attitude. We applied the FFT to computed gyro 
rates samples that were of different lengths (15 and 30 min) and different sampling times (1.152 and .128 sec, 
respectively). The units of the PSD are rad 2 /sec 2 . 

In Figures 5, 6, and 7, we have the power spectra density function from the quiescent period. The spectra indicates 
that the noise level of the PSD is about 10" 12 rad 2 /sec 2 . This implies a white noise level corresponds to a digital rate 
of about 0.2 arcsec/sec. On the X and Z-axis, there are four peaks that are in common to both plots, namely, .24, 
.97, 1.95, and 2.93 Hertz (Hz). The major difference between these axes is that the peak at .24 Hz on the X-axis is 
10 times larger than one on the Z-axis. The Y-axis, which is the pitch axis, is the axis about which UARS rotates 
once an orbit (1 rpo) and its power spectra does not have a well-defined peak at 0.24 HZ. The pitch axis does, 
however, have several frequencies of interest. The frequencies are as follows: 0.479 , 0.956 , 1.43, 1.95, 2.44, 
2.93., 3.43 Hz. 

Data from August 8, 1993, in which only the solar array is operating, shown in Figures 8, 9 , and 10, indicates a 
large increase in the amplitude on all three axes of at least three orders of magnitude at .24 and .956 Hz. In 
addition, there is a small peak at 2. 15 Hz in the power spectra of the X and Z axis. The pitch axis does not contain 
this frequency. 

Data from August 11, 1993, in which all the science instruments are turned on, does not indicate any new 
frequencies with a finite width. Likewise, the data from February 6, 1994, does not indicate any significant new 
frequency information other than a slight increase in the noise level. Figures 1 1 and 12 are representative samples 
of the data from the period when the science instruments were turned on and off from nominal operations. 

Discussion of Results 

The immediate goal of the power spectra analysis, which was to determine whether the method can be used to 
identify spacecraft motion due to the solar array and science instrument, has been obtained. The signature of the 
solar array consists of large amplitude noise at 0.24 and 0.967 Hz. Additionally, there is a small peak at 2. 15 Hz. 
There does not seem to be a noise signature associated with the operation of the science instruments. 

However, further analysis is needed to explain the presence of a large number of peaks with very narrow frequency 
bandwidth. The occurrence of these frequencies may be related to bandwidth. It is well known that high-frequency 
resolution (small bandwidth) leads to large variances of the estimate of power spectra while low-resolution (wide 
bandwidth) produces a stable estimate (Reference 7). On the other hand, presence of these frequencies may be 
related to how we compute the gyro rates. The basic time step in both these cases is .128 sec. Finally, the sharp 
peaks may indicate aliasing or over sampling The presence of increased noise levels near the Nyquist frequency is 
the standard indicator of aliasing 

A comparison between the UARS and EUVE power spectras was performed to determine if EUVE gyros had a 
similar signature to UARS. The calculation of EUVE gyro rates is completely different from that of UARS The 
data in Figure 13 indicate well-defined signal peaks at 0.976, 1.953, and 2.93 Hz. The EUVE power spectra in 
Figure 13 represents FFT of 8,192 points from a time series of 16,384 points Since the Telemetry Processors (TPs) 
for EUVE and UARS are so different, but mathematically equivalent, the weight of evidence indicates that the 
abovementioned peaks may be characteristic of the Teledyne DRIRU II gy roscopes and/or multimission spacecraft 

The presence of similar peaks in both power spectral densities rules out computational error in the way rates are 
computed but does not rule out aliasing. The FFT computer program does not calculate FFT above the Nyquist 
frequency of the time series, but the time series has not been filtered to remove frequencies above the Nyquist 
frequency. Note that every point in the time series contributes, in principle, to the amplitude of every harmonic. 
Further work needs to be done to construct a time series filter that would provide a definitive answer about the 
significance of the narrow peaks 
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Figure 5. Power Spectra of the X Axis: 920604: (Sample = .128 sec)(N - 8192) 
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Figure 6. Power Spectra Y Axis: 920604 (Sample time = -128 sec)(N = 8192) 
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Figure 7. Power Spectra Z Axis: 920604 (Sample time = .128 sec)(N = 8192) 
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Figure 11. Power Spectra X-axis: 930811 (Sample time = .128 sec) (N - 8192) 
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Figure 12. Power Spectra X-axis: 940206 (Sample time = .128 sec) (N = 8192) 
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Figure 13. EUVE Coefficients for Z Channel (7/2/92) (N = 16384) 
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